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We examine if there is a lower bound on the detection cross section, cr^ 



for the neutralino 



dark matter in the MSSM. If we impose the minimal supergravity boundary conditions as well as 
the "naturalness" condition, in particular mi/2 < 300 GeV, we show that there is a lower bound 
of CTx-p > 10"''® cm^. We also clarify the origin for the lower bound. Relaxing either of the 
cissumptions, however, can lead to much smaller cross sections. 



I. INTRODUCTION 

Supersymmetry is considered to be a compelling ex- 
tension to the Standard Model for several reasons. For 
example, it stabilizes scalar masses against radiative 
corrections, allowing theories with fundamental scalars 
to become natural. For a review see [Q . Supersymmetry 
also weighs in on the dark matter problem: stars and 
other luminous matter contribute a small fraction of the 
critical density, ^lum = (0.003 ± O.OQl)/i~^, while the 
amount of matter known to exist from its gravitational 
effects (both at galaxy and cluster of galaxies scales) is 
much larger, D,m — 0.35 ± 0.07 (Ref. M). Furthermore, 
most of the missing matter seems to be non-baryonic 
in nature. The experimental motivation behind the 
dark matter problem and different search strategies are 
discussed in more detail in |, |, §. 

In supersymmetric models, R-parity is often imposed 
to avoid weak-scale proton decay or lepton number 
violation. Imposing this symmetry also yields an ideal 
fermionic dark matter candidate. Namely, in supersym- 
metric models with R-parity, the lightest supersymmetric 
particle (LSP) is stable and it could conceivably make up 
a substantial part of the dark matter in the galactic halo. 

Here we investigate the direct detection of such a 
particle. There have been many such studies in the 
literature [|[ |, 0, |, |9[ |l0). More recently, there have 
been discussions of numerous variables that can effect 
direct detection. These studies include an investigation 
of the effect of the rotation of the galactic halo the 
effects of the uncertainty of the quark densities within the 
nuclei ||l2, 13, 14j, possible CP violation |l^,|T^, and non- 
universality of gaugino masses [|l^ . Here we attempt 
to address the following question: Is there a minimum 
cross section for the elastic scattering of neutralinos off of 
ordinary matter? Naively, it would seem that a judicious 
choice of parameters might allow a complete cancellation 
between different diagrams. After all, the parameter 
space is very large in the general Minimal Supersym- 



metric Standard Model (MSSM), and even for a very 
restrictive framework such as the minimal supergravity 
(mSUGRA), the number of parameters is still quite large. 
We will show that there nonetheless exists a minimum 
cross section in the mSUGRA framework. However, 
we will also show that this result strongly depends on 
the assumptions of the framework, such as unification 
of different parameters at the GUT scale, radiative 
electroweak symmetry breaking and naturalness. 

This argument is of great importance when considering 
the upcoming direct detection experiments. For the 
mSUGRA framework, one expects that the future am- 
bitious direct detection experiments can explore most of 
the parameter space. However, we find that the detection 
picture is not quite as rosy for a more general MSSM 
framework. 



II. DEFINITIONS AND APPROACH 

We adopt the following notation for the superpotential 
and soft supersymmetry breaking potential in the MSSM: 



W 



Vsoft 



-u*j^AuhuqiH^ - BfiHlH^ + h.c.) 



(1) 



-d*^Mj,dR 



--{MiBB + IvhWW + IvhgY)- 



(2) 



Here the Ks are Yukawa couplings, the A's are trilinear 
couplings, the Mq^u^d^l^e are the squark and slepton 
mass parameters, the Mi_2.3 are gaugino mass parameters 
and , w/Za i Mi and B are Higgs mass parameters. 
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The i and j are SU{2)l indices, and are made explicit, so 
as to make our sign conventions clear. SU (3) indices are 
suppressed. In the R-parity invariant MSSM the LSP is 
usually a neutralino - a mixture of bino, neutral wino and 
two neutral higgsinos. In our notation, the neutralino 
mass matrix reads 

(Ml ~mzse„,ci3 +mzsg„si3 

—rnzsg^^^cp ~\-irLzCe„Cis — 

\ -TnzSe„Sj3 -mzce^sp -fi 

■(3) 

Here — sin/3, = cos/3, sg„ = sinOw, and ce„, = 
cos 9w The physical states are obtained by diagonalizing 
this matrix. The lightest neutralino can be written in the 
form: 



Xi 



NnB + N12W3 + NisH^ + (4) 



We are interested in spin independent scattering of 
neutralinos off of ordinary matter. This contribution 
dominates in the case of detectors with large nuclei, 
such as Ge As discussed in the literature, in 

most situations the dominant contribution to the spin 
independent amplitude is the exchange of the two neutral 
Higgs bosons, although in some cases the contribution of 
the squark exchange and loop corrections are substantial. 
The relevant tree-level diagrams are shown in Figure |l|. 

We use the DarkSUSY package ^ in this study. This 
code has the following inputs: Mi^2,3, the ratio of 
the vevs of the two Higgs bosons (tan/3 = V2/V1), the 
mass of the axial Higgs boson (rriA), the soft masses of 
the sparticles {Mq^jj^d^l^e) and the diagonal components 
of the trilinear coupling matrices {Ae,d,u)- All inputs 
are to be supplied at the weak scale. DarkSUSY then 
calculates the particle spectrum, widths and couplings 
based on the input parameters. It evaluates the cross 
section for scattering of neutralinos off protons and 
neutrons, following Ref. [nOl. It also evaluates the relic 




FIG. 1: The leading diagrams for direct detection. Note that 
there is also a u-channel diagram for squark exchange. There 
are also diagrams where the neutralino scatters off of gluons 
in the nucleon through heavy squark loops. 



density of the neutralinos for the given input parameters 
following |2^, which includes the relativistic Boltzmann 
averaging, sub-threshold and resonant annihilation and 
coannihilation processes with charginos and neutralinos. 
Furthermore, DarkSUSY checks for the current con- 
straints obtained by experiments, including the 6^5 + 7 
constraint p2[ . 



III. MSUGRA FRAMEWORK 

A. Definition of the Framework 

In this section, we briefly outline the mSUGRA frame- 
work. In mSUGRA, one makes several assumptions: 

• There exists a Grand Unified Theory (GUT) at 
some high energy scale. Consequently, the gauge 
couplings unify at the GUT scale. The value of the 
couplings at the weak scale determines the GUT 
scale to be ~ 2 X 10^^ GeV. The gaugino mass 
parameters also unify to mi/2 at the GUT scale. 

• Other unification assumptions are: the scalar mass 
parameters unify to a value denoted by mo and the 
trilinear couplings unify to Aq at the GUT scale. 
Using the MSSM renormalization group equations 
(RGEs) we evaluate all parameters at the weak 
scale. We choose to do this using the one loop 
RGEs that can be found, for example, in |23] or 



• Radiative electroweak symmetry breaking 
(REWSB) is imposed: minimization of the one- 
loop Higgs effective potential at the appropriate 
scale fixes and tua (we follow the methods 
of Refs. 1^ and |^). For completeness, we 
reproduce the equation for /i^ at tree level: 



'Hi 



tan' 



/3 



1 



tan^ /3 - 1 



-Mi 



(5) 



With these assumptions, the mSUGRA framework 
allows four free parameters (mo, TOi/2, ^0 and tan/3). 
Also, the sign of fi remains undetermined. Starting 
with these parameters we determine all of the input 
parameters for the DarkSUSY code. We allow the free 
parameters to vary in the intervals 

< mi/2 < 300 GeV, 95 < ?7io < 1000 GeV, 
-3000 <Aa< 3000 GeV, 1.8 < tan/3 < 25. (6) 

The upper bounds on these parameters come from the 
naturalness assumption: one of the reasons for using 
supersymmetry is its ability to naturally relate high 
and low energy scales; as a result, no parameter in the 
theory should be very large. Moreover, the upper bound 
on mi/2 makes our data set insensitive to the stau- 
neutralino coannihilations, which are not included in the 
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calculation of relic density performed by DarkSUSY. We 
have explicitly checked that there are no models in our 
data set which would be cosmologically allowed only if 
the stau coannihilations were included. We will also see 
later that expanding the upper bound on TOi/2 has serious 
consequences regarding the lower bound on the elastic 
scattering cross section. Also, the low value of tan/3 is 
set by the requirement that the top Yukawa coupling does 
not blow up before the GUT scale is reached. 

Before we present the detailed analysis of cross section, 
a few remarks are in order. First, the b — > s+7 constraint 
eliminates large portions of the fi < parameter space, 
in agreement with |2^, |28|. Second, for both /j, > 
and /X < 0, we find no higgsino-like LSP models that are 
cosmologically important, in agreement with p8| . 

We plot the variation of the spin independent cross 
section versus the neutralino mass in Figure |^. Note 
that the complete allowed region is split into two parts 
by the annihilation channel into W~^W~ , which affects 
the relic density of neutralino. We consider two relic 
density constraints. Since the present observations favor 
the Hubble constant h = 0.7 ± 0.1 and the total matter 
density ^Im — 0.3 ± 0.1, of which baryons contribute 



0.236. 



0.02, 



we consider the range 0.052 < fly-t 



< 



However, we also examine effects of relaxing the 

h"^ < 1. The upper 



relic density constraint to 0.025 < H.-^ 
bound on cr^-p (of the theoretically allowed regions) 
comes from the lower bound on the relic density. The 
lower bound on comes from the existing constraints 
from the accelerator experiments. The upper bound 
on is a combination of the upper bound on relic 
density and of the bounds on the free parameters. The 
lower bound on cr^_p is not yet well understood, and 
it is the subject of this paper. Figure § also includes 
some recent and future direct detection experimental 

33|| . Note that the parameter 
) corresponds to the region of 
plane bounded by the two closed solid 



results |3l|, |3 

space defined by Eq. 
a 



cr^-p > 10 cm for these models, 



'x-p 

lines. Therefore 
or equivalently, assuming a ^"^Ge target, the dark matter 
density pD — 0.3 GeVc^^cm^^, the WIMP characteristic 
velocity vq — 230 km s^^ and following Ref. Q, the 
event rate R > 0.1 ton~^day~^. Hence, the most 
ambitious future direct detection experiments may be 
able to explore a large portion, if not all, of the these 
models. 




WIMP Mass (GeV) 

FIG. 2: The cross section for spin-independent x-proton 
scattering is shown. Current accelerator bounds, including 
6 — > s + 7 are imposed through the DarkSUSY code. The 
top dark region is the DAM A allowed region at 3cr CL. The 
dashed curve is the DAMA 90% CL exclusion limit from 1996 
(obtained using pulse-shape analysis) . The thinner solid curve 
is the current CDMS 90% CL exclusion limit, the thicker solid 
curve is the projected exclusion limit for CDMS II experiment 
and the dotted curve is the projected exclusion limit for 
the GENIUS experiment. The light shaded regions denote 
mSUGRA models passing the 0.025 < flh^ < 1 constraint 
and with the upper bound mi/2 < 1 TeV. As discussed in the 
text, this region shows that including stau coannihilations 
into the relic density calculation has a dramatic effect on 
the lower bound on cr^-p for mi/2 > 300 GeV {M^ > 120 
GeV). Restricting the relic density further to 0.052 < Qh^ < 
0.236 yields the darker regions within the lighter regions and 
restricting also the upper bound on mi/2 < 300 GeV (Eq. 
(^)) gives the regions bounded by the closed solid curves. 



therefore, concentrate on the Higgs boson exchange and 
will postpone the discussion of squark exchange to the 
end of this section. 

The contribution of Higgs boson exchange can be found 
in the literature [| |3[ H || |6[ |7|. It is of the 
following form: 



B. Results and Analysis 



<Jh,H -\{.fu+fc + ft)A" + ifd + fs + .MA^'l 



d\2 



(7) 



As mentioned above, the dominant contribution to 
spin independent elastic scattering is usually the Higgs 
boson exchange. Figure ^ illustrates this relationship 
within our results. The nearly perfect 45 degree line in 
the figure indicates good agreement between the total 
cross section as evaluated by DarkSUSY and the cross 
section calculated including the exchange of Higgs bosons 
only (in the approximation explained below). We will. 



where /„ « 0.023, « 0.034,/, « 0.14,/, = /( = /,« 
0.0595 parametrize the quark-nucleon matrix elements 
and 



A" 



92 



Fii cos an Fh sin an 
AMw V sin f3 mj^ sin /3 

ffl / Fh sin an _^ Fh cos an ' 
AMw ^ mf^ cos (3 mjj cos (3 . 



(8) 



(9) 
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0.025< Q < 1 

Accelerator bounds not applied 




H,h exchange only 



10 

(cm') 
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FIG. 3: The cross section for spin-independent x-proton 
scattering: the complete (DarkSUSY) calculation is shown 
on the y axis and the contribution to the cross section from 
the Higgs bosons exchange alone is shown on the x axis. A 
relatively conservative relic density cut is applied, 0.025 < 
n/i^ < 1. 




- Case 3 models 
_2{) - accelerator bounds not applied 
0.025 < £2 h'^ < 1 



-1000 -800 -600 -400 -200 200 400 600 800 1000 
H (GeV) 

FIG. 4: This figure illustrates the regions of Mi — /i plane 
(as defined by the Eq. (^)) in which the different cases 
discussed in the text can be satisfied. For the models shown, 
the accelerator constraints were not applied, but a relatively 
conservative relic density cut (0.025 < Q,h^ < 1) was applied. 
Note that almost all of the models of Case 3 shown in this 
plot are ruled out by the accelerator constraints. Also, there 
are no models satisfying the Case 5 for the parameter space 
defined by Eq. (H) . 



and 

Fh = {Ni2 - A^ii ta.ii9w){Nii cos an + A^is sinan) 
Fh = {Ni2 - Nil tan0w)(iVi4 sin an - N13 cosai^|0) 

The A^'s are the coefHcients appearing in Eq. (jj) and an 
is the Higgs boson mixing angle (defined after radiative 
corrections have been included in the Higgs mass matrix) . 
Note that there is an upper bound on the light Higgs 
boson mass in the MSSM, given by mh < 130 GeV [ p8[ . 
A" represents the amplitude for scattering oS^ an up-type 
quark in a nucleon, while A"^ represents the amplitude for 
scattering off a down-type quark in the nucleon. Since all 
models we generate have a bino-like neutralino, /i > Mi 
and fi > Mz- Then, following Ref. we can expand 
the A^ii's out in powers of We reproduce their result 
here: 



Nil 
N12 

Ni3 
Ni4 



1, 



sin 26* 



w 



M, 



2 ^ (1 - Mi7/x2) M2 - Afi 
Mz 1 



sin 2/3 - 



(11) 
Ml' 



M l-Mf/f,^ 
Mz 1 



sin dw sin (3 
sin dw cos P 



1 



M- 



- cot (3 



1 



M: 



- tan/3 



In this approximation, we can identify five possible sit- 
uations in which the neutralino-proton elastic scattering 



cross section is relatively small. Figure Q illustrates which 
regions of the Mi-fi plane satisfy the conditions of the five 
cases. In the following, we will use the approximation of 

Eq. dnl). 

• Case 1: iVi2 — Nn ta.n6w = would make both Fh 
and Fh vanish. Intuitively, this is reasonable since 
this condition implies that the neutralino is a pure 
photino, and the tree level Higgs coupling to the 
photino vanishes. Using Eq. (|l l|) , we may rewrite 
this condition as 



H{M2 - Ml) = -Ml cos^ Ow ( sin 2/3 



Ml 



(12) 



Since M2 > Mi > 0, the last equation can be 
satisfied only if /i < 0. If we use the GUT 
relationship betweeen Mi and M2, we get 



f,Mi ('^£^!!^ _ 1) = -Mlcos^e^ ('sin2/3+ ^' 
5 sin 9w J \ M / 

(13) 



or, equivalently, for and Mi in units of GeV 
-6449 sin 2/3 



Ml 



6449 



(14) 



As sin 2/3 ranges from to 1, Eq. (Q) spans the 
dashed regions marked 'Case V in Figure 0. For 
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sin 2/3 = 1, Eq. @ implies « Mi < 40 GeV. 
This can be read directly from Figure ^. Note that 
such low values of are excluded by the relic 
density constraint and the current experimental 
limits (as shown in Figure ^ and Ref. The 
constraint on M-^ is even stronger for other values 
of sin 2/3, so we conclude that this condition cannot 
be satisfied in the mSUGRA framework. 

Case 2: Nij^snvan — Ni^cosan = would make 
only Fh vanish. In our approximation, this condi- 
tion translates into 



Ml I ^ = — cot an — cot /?. 



(15) 



To understand the meaning of this condition better, 
we will use the tree level relationship between an 
and /?: 



]\'f^ 

cot2aH = fc cot 2/3, k = f. 



(16) 



which, after some trigonometric manipulations, 
yields 



k f I 

cot an — — I ■; — tan p 



2 Vtan/5 



fk^/ 1 
y T Vtair 



+ tan2/3~2) +1. (17) 



Since both terms on the right hand side of Eq. (|l7| ) 
are negative because tan/? > 1, the minimum of 
I cot a// 1 — 1 occurs when fc = (or, equivalently, 
when niA = Mz). Then, since tan/3 > 1.8, the 
condition of Eq. ( [l5| ) becomes 



Ml 



> 0.5. 



(18) 



At this point, we would like to formulate a re- 
lationship between fj, and Mi resulting from the 
RGEs and REWSB assumptions. In Ref. [||, 
an approximate solution (based on the expansion 
around the infrared fixed point) to the RGEs for 
the Higgs mass parameters, ttihi ^^'^ '^H2i is 
presented. Assuming that the value of the top 
Yukawa coupling is relatively close to the infrared 
fixed point, we can write: 



riQ + 0.5toi 
-0.5too ~ 3.5r7ii 



/2' 



/2- 



(19) 



These equations, coupled with Eq. (^, yield a value 
for fi^ in terms of toi/2, wq, and tan/3. Using the 
GUT relationship: 



Ml = mi/2 



ai{rnz) 

OiGUT 



(20) 



we get a roughly linear relationship between /i and 
Mi: 



Ml = (0.3|/i| - 60) ±40. 



(21) 



The spread ±40 comes from the variation in mg 
and tan /3 and the linearity breaks down somewhat 
at the low values of Mi . The empirical relationship 
that we obtain from running the code (see Figure 
Wj is very similar: 



Ml = (0.3|Ai| -40) ±25. 



(22) 



The smaller spread comes from the application of 
the relic density cut and the current experimental 
limits. In any case, we conclude that Mi/fi > 0.3 
is not allowed in the mSUGRA framework. This is 
in conflict with Eq. ([l^), implying that the Case 
2 cannot be satisfied. Note that the tree level 
relationship in Eq. (^6|) is altered at higher orders, 
but we have checked that this does not affect the 
final conclusion. 

• Case 3: Nucosan + Ni^smaH = would make 
Fh vanish. Manipulation of this condition using 
Eq. (|ll|) yields 



Ml/ ^ — tana// — cot/3. 



(23) 



Figure ^ shows that when this condition is (approx- 
imately) satisfied, the contribution of the heavy 
Higgs boson exchange to the elastic scattering cross 
section is the largest. Hence, even if this condition 
is satisfied, <Jx-p cannot be arbitrarily small due 
to heavy Higgs boson exchange (along with other 
channels such as the squark exchange) . Of course, 
this assumes that there is some upper bound on 
the heavy Higgs mass, which is true simply because 
the parameter space is bounded. Note, also, that 
since ianan < 0, the condition for the vanishing of 
the light Higgs boson contribution can be satisfied 
only for fx < 0. It is interesting, however, that most 
of the models in our data set that approximately 
satisfy this condition are excluded by the bounds 
on Higgs boson masses. 

• Case 4: There is one more way of making both F^ 
and Fh small, and that is by making fj, very large 
{Ni2, Ni3, and A^i4 all contain /i in denominator). 
However, this possibility is limited by naturalness 
assumption: fi is kept below « 850 GeV by the 
upper bound we have chosen on mi/2- Hence, in 
mSUGRA framework, the naturalness assumption 
also keeps the cross section from vanishing. 

• Case 5: We now consider the possibility of complete 
cancellation of terms in the Equation (0). Let 
us examine the relative signs of F^, Fh and fi. 
Since tan/3 > 1.8 and Mi/ji k, 0.3, it follows from 
Eq. (|ri|) that |A^i3| > |/Vi4|. Then, since cot a// < 
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FIG. 5: The proton-neutralino scattering cross section for 
models where the condition in Case 3 is approximately satis- 
fied (in particular, models satisfying 0.8 < tan a^'^-cot < 1/2 
are denoted by o's). The complete (DarkSUSY) calculation is 
shown on the y-axis and the contribution to the cross section 
due only to the exchange of the heavy Higgs boson, is shown 
on the a;-axis. A relatively conservative relic density cut is 
applied, 0.025 < fih^ < 1. Most of the models denoted by 
o's in this plot are excluded by the accelerator bounds on the 
Higgs boson masses. 

— 1 (and sin an < 0), the A^i3 term dominates 
over the A^i4 term in Fh (Eq. ^)). Hence, 
Fh / 1^ > always, consistent with the analysis of 
Case 2 above. The situation is somewhat more 
complicated in the case of Fh- For > 0, following 
a similar analysis we get Fh/fi > 0. Then, the 
interference between the two terms in A" (Eq. (||)) 
is destructive and the interference between the two 
terms in A''- (Eq. (||)) is constructive. Furthermore, 
in Eq. (|^) we see that A" gets multiplied by a 
much smaller form factor than A*^. As a result, A"^ 
strongly dominates over A" in Eq. (^, preventing 
<7x-p from vanishing. 

On the other hand, if /i < and ii fi > —Mi tan/?, 
iVi4 can change sign relative to fi. This change 
of sign can propagate through Equations (^0|), (0), 
and (0), so that the A" and A"^ terms in the Eq. ^) 
are of opposite sign. If the parameters are tuned 
properly, a complete cancellation of these terms can 
be obtained. However, this cancellation happens 
only if Mx(w Mi) > 120 GeV, which is not allowed 
due to the upper bound on toi/2. We conclude 
that in the parameter space defined by Eq. (||) , A''' 
always dominates over A", so cr^-p does not vanish. 



If the upper bound on mi/2 is relaxed to 1 TeV, 
the situation changes significantly. The complete 
cancellation discussed above is now allowed to 
happen. As shown in Figure ^, the lowest values 
of cr;^_p are reached exactly when the A'^ and A" 
terms cancel each other. As shown in Figure H the 
stringent relic rensity cut 0.052 < ilh^ < 0.236 
rules out all of these models. However, we do 
not include the stau-neutralino coannihilations in 
the calculation of the relic density, so this result 
should be taken with caution (see, for example, 
g 0). In fact, our data set contains models 
in which the stau-neutralino coannihilation could 
make a difference. For this reason, the Figure ^ 
also contains the allowed region for a conservative 
relic density cut 0.025 < flh^ < 1. Clearly, this cut 
allows models with very low cr^-p values. 
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FIG. 6: The absolute values of the and the A" terms 
of Eq. (0) are shown. Only < ^ (for which M < 0) 
are shown, satisfying a conservative relic density constraint 
0.025 < Jl/i'^ < 1, passing the accelerator bounds, and with 
"11/2 < 1 TeV. The o's denote the lowest CTx-p models (< 
10~** cm^). Clearly, the lowest values of a-^-^ are achieved 
when the A^ and A^ terms cancel each other out. 

With the above discussion in hand, let us go back and 
consider the squark exchange. The complete calculation 
of the squark exchange contribution is fairly complex. 
However, good insights can be gained by making several 
simplifying assumptions. First of all, the contribution of 
the squark exchange can be roughly approximated by the 
contribution of the exchange of the u, d, and s squarks. 
In this case, the contribution of the squark exchange to 
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the cross-section can be written as 

<^q ^ \ fuBu + fdBd + fsBs 



(24) 



where fi is as defined above, and the Bj represent the 
amplitude for scattering off of a quark of type j in the 
nucleon. Furthermore, in the following considerations 
we neglect the left-right mixing in these light squarks. 
This should be true over a large class of models, as 
the off-diagonal elements in the squark mass matrix 
are proportional to the corresponding quark mass. Let 
us also neglect the mass splitting of the two different 
squarks. Also, since fs 3> /u, fd, we can neglect all but 
the Bg term. In this approximation, following Ref. p5[ , 
we can write: 



Bs = — 



1 



1 



4m, Mi 



[2C1C2 - 2C1C3 



(25) 



where we have defined the following: 



C2 
C3 



g2msNi3 
2Mw cos 13' 

eQyi 
eQyi 



cos 6w 

92 

cos 6w 



2/2 [73 - Qsin^ 6w], 
2/2(9 sin^ 9w- 



(26) 



Note that Ci represents the coupling of the down type 
quark to the Higgsino portion of the ncutralino. C2 
and C3 represent the couplings of bino to the left and 
right handed quark, respectively. Here, T3 is the SU{2) 
quantum number of the squark in question, Q is the 
charge, yi denotes the photino fraction of the neutralino, 
while 1/2 denotes the zino fraction. They are given by: 



2/1 = A^ii cos6'w + A^i2 sin6'vi/, 
2/2 = - A^ii sin 9w + N12 cos 9 



(27) 



After approximating 2/2 ~ — sin 9w and using 
tan^n/ = (/'/(/, a brief and straight-forward calculation 
yields a simple expression for the amplitude due to the 
exchange of the strange squarks: 



-929 



N- 



13 



8Mw cos P M? - Ml ■ 



(28) 



Furthermore, we can write the masses Ms and in 
terms of the input parameters of mSUGRA. This is 
because the Yukawa couplings can be neglected in the 
RGEs. Following the methods described in Ref. and 
using the Eq. (|2(]| ) we can write 



Mi 



Ml 



5.8to? 



/2- 



(29) 



Using Eqs. ^ and ([2^), we can compare the squark 
exchange to the light Higgs boson exchange: 



j^d 2 sina//(A^i4 cosa/f -|- A^i3 sina/f)(mQ 



Here we have only kept the contribution of the strange 
quark to the Higgs exchange amplitude [A'^) as well. 
Note that in general, the light Higgs boson contribution 
will dominate. As expected, this is basically due to the 
fact that squarks are in general heavier than the lightest 
Higgs boson. The squark exchange can be important only 
if the contribution from the exchange of the Higgs bosons 
is fine-tuned to be very small. 



IV. GENERAL MSSM FRAMEWORK 

A. Definition of the Framework 

In this framework we relax our assumptions. We 
keep the unification of the gaugino masses, but we 
drop the requirements that the scalar masses and the 
trilinear scalar couplings unify. In addition, we drop the 
REWSB requirement (i.e., we take m^^^ ^ as independent 
parameters from mo). We assume that all scalar mass 
parameters at the weak scale are equal: msq. This 
assumption is made in order to simplify the calculation, 
and it should not affect the general flavor of our results. 
Of all trilinear couplings, we keep only At and A^ and 
we set all others to zero. Then, the free parameters 
are fi, M2, tan f3,mA,'rnsq, At, Ai,. We also relax the 
naturalness assumption, allowing the free parameters to 
have very large values. Besides the relatively uniform 
scans of the parameter space, we also performed special 
scans in order to investigate the different conditions 
mentioned in the previous section. The free parameter 
space is then: 



- 300 TeV < M < 300 TcV, 
95GeV <mA< 10 TcV, 

-3 < ^ < 3, 
rrisq 



< M2 < 300 TeV, 



200 GcV < 



^sq 



< 50 TeV, 



1.8 < tan /3 < 100. (31) 



Again, a few comments are in order. First, in this 
framework we observe higgsino-like (as well as bino-like) 
lightest neutralino. In agreement with the Ref. p2[ , 
we find very few light higgsino-like models, which will 
probably be explored soon by accelerator experiments. 
Most of the higgsino-like models (with gaugino content 
Zg < 0.01) have M^^ > 450 GeV, implying very large 
values of mi/2. In particular, in higgsino like models 
Ml > « M^; our resuhs give Mi > 700 GeV or, 
equivalently, mi/2 > 1700 GeV, which can be considered 
unnatural. For these reasons, we choose not to analyze 
the higgsino case. Second, 6 — > s -I- 7 is less constraining, 
but our results are still consistent with Refs. j2^, psj . 

We present the plot of a^-p versus M^ in this frame- 
work (Figure 0) . We do not pretend that Figure |^ reflects 
all points accessible in a general MSSM. However, it does 
serve to show some generic differences from the mSUGRA 
case. Namely, we can obtain much larger values for the 
neutralino mass because of the size of the parameter 
(30) space. In addition, the lower bound on c^^p is also much 
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FIG. 7: The cross section for spin-independent x-proton 
scattering in the general MSSM framework is shown. A 
relatively conservative relic density cut is applied, 0.025 < 
^ih^ < 1. Note also that the constraint on the gaugino 
fraction, Zg > 10, is applied. Current accelerator bounds, 
including 6 — > s-l-7 are imposed through the DarkSUSY code. 



lower than in the mSUGRA case. We discuss the specifics 
of this below. 



B. Results and Analysis 

We concentrate only on the bino-like lightest neu- 
tralino. All results in this section are presented with 
this assumption in mind. In particular, we demand 
Zg = {N^^ + Nl^)/[Nl^ + Nl^) > 10. In this case, we can 
rely on the same approximations we used in the previous 
section. In particular, the expansion of Eq. ( pd] ) is valid. 
So, we can revisit the 5 different cases explored in the 
previous section. 



• Case 1: N12 — Nu ta,n9w = cannot be satisfied 
because, as in the mSUGRA case, it implies < 
40 GeV, which is ruled out by the relic density 
cut and the experimental limits (as shown on 
Figure 0). The models that get close to satisfying 
this condition have very low contributions due to 
the Higgs bosons exchange, so this is one of the rare 
situations where the squark exchange is important. 

• Case 2: N14 sin an — cos an — Ois now possible 
to satisfy because /i and Mi are not related. 
We indeed observe that the heavy Higgs boson 
exchange contribution is very small in this case, but 
since the light Higgs boson exchange dominates, 
Ux-p kept relatively high in value. 

• Case 3: A^i4 cos an + Ni^ sin an — Ois also possible 
to satisfy. As in the mSUGRA case, the light Higgs 



boson exchange is small and heavy Higgs boson 
exchange dominates. Unlike in the mSUGRA 
case, the accelerator bounds do not rule out these 
models. 

• Case 4: Since the naturalness constraint has been 
relaxed, is allowed to have very large values. 
Then N12, and Nu can be driven small, which 
in turn would make the Higgs boson exchange 
contribution small. Intuitively, large implies 
that the neutralino is a very pure bino, for which 
the Higgs boson scattering channels vanish. If the 
squark masses are kept large as well, the squark 
contribution will be small too, making the total 
elastic scattering cross section very small. This is 
illustrated in Figure |^ - the lowest values of cr-^-p 
are obtained for the largest values of \fi\. 

• Case 5: As discussed in the mSUGRA case, if 
(■^tan/3) is negative and sufficiently large, A^i4 
can change sign. Following through Eqs. ( [To| ) and 
(|), this effect can induce destructive interference 
between the and A'^ terms in Eq. (0) and 
cause the overall cr^-p to vanish. We observed this 
cancellation in the mSUGRA case for > 120 
GeV, and we also observe it in the general MSSM 
case. In Table 1 we present some of the models in 
which this kind of cancellation takes place. Besides 
the models at large M^, we also observe models 
with relatively low Af^(Ri Mi) and very low values 
of (Jx-p (see rows 3 and 4 of Table 1). These models 
were not allowed in the mSUGRA framework due to 
the Ml — fi relationship (determined by the RGE's 
and the REWSB assumption), which does not exist 
in the general MSSM framework. 

In other words, in the general MSSM framework it 
is possible to obtain low values of cr^-p several 
different ways: either by tuning parameters to suppress 
the contribution of the heavy or the light Higgs boson 
exchange, or by allowing parameters (such as /i) to be 
very large, which suppreses both the heavy and the 
light Higgs boson exchange channels, or by fine-tunning 
parameters to achieve a complete cancellation of terms 
in the Eq. (^. As a result, the lower bound on cr^-p 
vanishes and it is beyond reach of the present and the 
proposed direct detection experiments. 



V. CONCLUSION 

We summarize our results as follows. The main 
contributions to the cross section for spin independent 
elastic scattering of neutralinos off nucleons come from 
the exchange of the Higgs bosons and squarks. The 
contribution of the squark exchange is usually much 
smaller, being of importance only when Higgs boson 
exchange contribution is very small. 
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^ (GeV) 


Ml (GeV) 


msq (GeV) 


ruA (GeV) 


tan/3 


At/msq 


Ab/msq 


(Tx-p (cm^) 


-1794 


502 


3792 


1004 


10.1 


1.2 


2.5 


7.9 xlO--'^ 


-2109 


534 


3211 


1087 


11.8 


-1.3 


1.0 


8.4x10"^^ 


-195 


55 


2995 


1120 


10.2 


-2.0 


2.3 


2.1x10^^^° 


-182 


61 


2891 


1099 


7.0 


-0.6 


-0.1 


1.6x10"^° 


-274 


163 


325 


1944 


3.8 


0.6 


2.5 


7.8x10"^° 



TABLE I: Some of the models in the general MSSM framework with very low values of 
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FIG. 8: The cross section for spin-independent x-proton 
scattering in the general MSSM framework is shown. A 
relatively conservative relic density cut is applied, 0.025 < 
Q,h? < 1. Note also that the constraint on the gaugino 
fraction, Zg > 10, is apphed. Current accelerator bounds, 
including b — > s + 7 are imposed through the DarkSUSY code. 

We investigate different conditions which could lead to 
small Higgs boson exchange contribution. We find that 
in mSUGRA framework, with the free parameter ranges 
defined in Eq. (^), these conditions are not satisfied 
due to the relationship between parameters Mi and /i 
(coming from the unification and radiative electroweak 
symmetry breaking assumptions), the naturalness as- 
sumption (which keeps different parameters from be- 
coming very large) and the accelerator constraints. We 
find that the light Higgs boson exchange dominates over 
the other channels and it leads to cr^-p 10^"*^ cm^. 
Equivalently, this yields an event rate > 0.1 ton^^day^^ 
in ^"^Gc target, which could be within reach of the future 
direct detection experiments. However, if we expand the 



mSUGRA framework to allow larger neutralino masses, 
the lower bound on <Jx-p vanishes for those large masses 
due to the occasional complete cancellation of terms in 
the Eq. (0). 

In the more general MSSM models (as defined in and 
above Eq. (|3ll)), the situation is significantly different. 
The light and/or heavy Higgs boson exchange channel 
can now be suppressed either by tuning parameters to 
satisfy Cases 2) or 3), or by allowing the free parameters 
(such as to be very large. Moreover, the complete 
cancellation of terms in the Eq. ^ is now possible even 
at low values of because the relationship between 
Ml and ^ that existed in the mSUGRA framework due 
to radiative electroweak symmetry breaking, is relaxed 
in the general MSSM framework. As a result, the lower 
bound on cr^-p is much lower than in the mSUGRA 
framework and it is beyond reach of the current or 
proposed direct detection experiments. 
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